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We consider unseeded Type I second-harmonic generation in quasi-phase-matched (QPM) quadratic nonlinear ma- 
terials and derive an accurate analytical expression for the evolution of the average intensity. The intensity-dependent 
nonlinear phase mismatch due to the QPM induced cubic nonlinearity is found. The equivalent formula for the intensity 
for maximum conversion, the crossing of which changes the nonlinear phase-shift of the fundamental over a period 
abruptly by it, corrects earlier estimates by more than a factor of 5. We find the crystal lengths necessary to obtain an 
optimal flat phase versus intensity response on either side of this separatrix intensity. 



Since the observation of nonlinear phase-shifts in £x.- 
cess of 7T through cascading close to phase-matchingErtl, 
quadratic nonlinear or x^ materials have been of sig- 
nificant interest in photonicst! With the maturing of 
the quasi-phase-matching (QPM) techniqueo, in partic- 
eld poling of ferro-electric materials, 
and by quantum-well disordering in 
, the number of possible applications of 



ular by electric- 
such as LiNbOp 
semiconductors^ 



cascading in x^ materials has increased even more. 

An effect of QPM gratings is to generate cubic non- 
linearities in the equations for the average field due to 
non-phase-matched coupling between modesQ. This cu- 
bic nonlinearity. appears in QPM with linear and/or non- 
linear gratingsa, it can be focusing or defocusing, de- 
pending on the sign of the phase mismatch^, and its 
strength can be significantly increased by modulating the 
gratingp. Simulations of QPM systems have confirmed 
the presence of an intensity dependent nonlinear phase 
mismatchEij and of soliton propertiesola, that can only be 
described by including the cubic terms. The Kerr nonlin- 
earity is also known to distort the second-harmonic (SH) 
spectrum and introduce an intensity dependent nonlin- 
ear phase misipatch, as found both theoreticallyM and 
experimentally!! 2 ]. Thus, although the physical origin 
of the Kerr and the induced cubic nonlinearities is dis- 
tinctly different, they have qualitatively similar effects on 
CW waves. 

In this letter we focus on the intensity dependent non- 
linear phase-mismatch, which implies a finite separatrix 
intensity that can be used in efficient all-optical switch- 
ing, as, was shown for the Kerr nonlinearity in japled 
fibresEj and the QPM induced cubic nonlincarityEj. A 
distinct feature of the averaged QPM model with induced 
cubic nonlinearities is that there is not a one-to-one cor- 
respondance between the physical field and the averaged 
field. The separatrix intensity obtained with no SH seed- 
ing in the averaged model does therefore not accurately 
predict the real physical separatrix intensity of unseeded 



second-harmonic generation (SHG). We address this dis- 
crepancy and find the exact separatrix intensity for un- 
seeded Type I SHG in QPM samples. We further show 
that the averaged model gives also quantitatively inac- 
curate results for no seeding of the average SH, since a 
basic assumption in the averaging procedure is violated, 
and we find the optimum crystal lengths for using the 
induced separatrix intensity for all-optical switching. 

We consider a linearly polarized electric field E = 
e[E\{z) exp(i/ciz — iut) +£2(2) exp(itt2Z ~ i2cut) +c.c.]/2, 
propagating in a lossless QPM x^ medium under con- 
ditions for type I SHG. The dynamical .eauations for the 
slowly varying envelopes take the formtil 



idEi/dz - 
idE2 1 dz 



G(z) X iE*E 2 e lAk 
■G(z) X 2Efe- iAk * 



= 0, 



(1) 
(2) 



where E\(z) is the fundamental wave (FW) with fre- 
quency u> and wavevector fci, £2(2) is the SH with 
wavevector k%, Afc=fc2— 2fci is the wavevector mismatch, 
and Xj =tj Jdeg / (rijc) , with nj=n(juj) being the refractive 
index and d e s=X being in MKS units. The to- 
tal intensity I — i^o^il-Eil 2 + n 2 |£ , 2| 2 ) is conserved, 
where r/o— v^o/Vo is the specific admittance of vac- 
uum. The x^ susceptibility is modulated by the grat- 
ing function G(z) with unit amplitude and Fourier se- 
ries G(z)=cr^ n g„e mKZ , where g n =0 for n even and 
g n =2/ (inn) for n odd. This zero-average square- wave 
modulation is typical for QPM by domain inversion 
in ferro-electric materials, such as LiNb03. We define 
cr=sign(/i) so that G(z) is positive in the first domain. 

We further consider first order forward QPM with a 
short coherence length Ld^L c ^L, where L is the crystal 
length, Ld=ir/\n\ is the domain length, and L c —Tr/\Ak\ 
is the coherence length. Expanding the fields in Fourier 
series in the grating wavenumber k, 



£1 



w n (z)e mKZ , E 2 = a^n{z)e 



i(riK,—(3)z 



(3) 
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where the harmonics are small compared to the 
component, a simple first-order perturbation theory 
gives the dynamical equations for the slowly varying (on 
the scale of the domain length) average field 

idwo/dz + ipiw^vo + (71M 2 - j2\wo\ 2 )w = 0, (4) 
idvo/dz + flvo — ipiw\ + 2j2\wq\ 2 vq = 0, (5) 

where (3=Ak— k<k is the residual mismatch, pj=Xj2/ir, 
and 7j=x.j'Xi(l~8/7r 2 )/ft. Since L d ~L c <^L the sign of k 
is the sign of the mismatch, sign(«)=sign(Afe). Thus the 
QPM induced cubic nonlinearity can be both focusing 
and defocusing, depending on the sign of the mismatch 
Ak, just as the effective cubic nonlinearity obtained in 
the cascading limitta. The harmonics are given by 



UKWn^o = Xl9n-lW V , nKVn-iQ = X29n+lW , 



(6) 



The average model (0jq) conserves the total average 
intensity Iq = ^r?o(^i l^or + "2 |^o | 2 ) an d the quantity 
e = (1 — u)y^tsin(0 2 — 2<fii) + Cu + Du 2 . Here 4>\(z) 
and 4>2(z) are the phases of wq(z) and vq(z), respec- 
tively, and u(z) = ^ 170 ^2 1 (2;) | 2 //o is the fraction of 
average intensity Iq in the average SH. The parameters 
are D — §7ivJ/pi and C — — \(3/(pi\/~J) — |D, where 
^0 = ^Vo n 2J- In the physical regime, where the polyno- 
mial f{u) = -D 2 u 4 + (1 - 2CD)u 3 - (2 + C 2 - 2eD)u 2 + 
(l+2eC)u— e 2 has four real roots uq < u\ < 1 < U2 < U3, 
we solve Eqs. by quadrature and find the average 

SH intensity 

u(z) = [it3 sn 2 (rz|m) + uop]/[sn 2 (rz\m) +p], (7) 



where p = (113 
and r = ||7i| J^/Juz 



U\)/{u \ — Uq), nip — (U3-U2)/(U2-U0), 

u o){u3 — Mi) are real positive pa- 
rameters. The Jacobian Elliptic sn(z|m) function is peri- 
odic with period 4K(m), whepe.K(m) is the complete el- 
liptic integral of the first kindEj. The solution (|^) is thus 
periodic with the period 2K (m)/r and becomes aperiodic 
when u\=U2=\ (m=l), in which case 100% transfer of 
power to the SH is predicted by the average model. From 
/(1)=0 we obtain the separatrix C + D=e. 

The solution (g) was found earlier fqs-.no seeding of 
the average SH, i.e. with u(0)=uo=e=(J 1 Si. In this case 
I=Iq to lowest order and the separatrix C+D=Q corre- 
sponds to the physical intensity L=I® =\ r no n 2 >/.?, where 
Jg=— /3/7x- It was further prover£3 that when the sepa- 
ratrix is crossed the phase-shift of the FW over a period 
changed by exactly tt, a result that is also valid for our 
general solution ([?]). However, it was never investigated 
whether 1®=— ^o n 2/3/7i was an accurate prediction of 
the separatrix intensity in the physical system ([l||^) . 

As a typical sample we consider bulk LiNbC>3, a 
FW wavelength of 1.064/im, c? c ff =30pm/V, ni=2.2, and 
ri2=2.23, which gives the coherence length L c =8.9/im. 
In Fig. |] we show the nonlinear phase-shift of the 
FW versus the physical intensity / as predicted by 
the unseeded average model (||-||) with it(0)=i>o(0)=0 



and found by numerical integration of the correspond- 
ing physical system d3-^|) with initial condition de- 
termined from Eqs. (g) and (|^). With the domain 
length Ld—S.8pm the predicted separatrix intensity 
7 s °=231GW/cm 2 is 26GW/cm 2 too high compared to 
what is actually found in the physical system (Q-||) . This 
inaccuracy has not been observed before and is mainly 
due to the average SH being initially zero, which violates 
the requirement that the harmonics are small compared 
to the average field. 
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Fig. 1. Nonlinear phase-shift of the FW versus intensity 
I for bulk LiNbOs, as predicted by the averaged model (4-|) 
(dashed) and found numerically from the physical Eqs. (1-E) 
(solid). The QPM domain length is Ld=8.Sp,m and the crys- 
tal length is L=351um (174/im) for unseeded (seeded) SHG. 
The vertical dotted lines mark the separatrices 7 S and 1°. 



At this point we stress that no seeding in the average 
model implies a seeding of the physical SH, as given 
by the transformations (||) and (||), i.e. «o(0)=0 gives 
Ei(0)=w {0) and E 2 (0)=-i(o-p 2 / n)w 2 (0). The experi- 
mentally relevant setup is unseeded SHG with i?2(0)=0, 
for which Fig. [l] shows that the physical phase versus in- 
tensity curve is indistinguishable from the prediction of 
the average model (vq is now seeded and thus no assump- 
tions are violated). The physical separatrix intensity for 
unseeded SHG is 44GW/cm 2 , i.e., 5.3 times lower than 
the prediction /,? of the unseeded averaged model. Note 
that the averaged model without the cubic terms would 
not predict any separatrix. 

From Eqs. (§) and (|) with E 2 {0)=0 we find that 
vo(0)=i(p2/ K)3r where x—wa(0) is real and Ei(0) — 
x+ (pip 2 /K 2 )x 3 . Then e—n 2 p2\ /r J / \nin) to lowest order 
and thus J s =—(1—8/tt 2 )(3/ji. This gives the physical 
separatrix intensity for unseeded SHG 
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(8) 



since again I=Iq to lowest order. For Ld— 8.8/im this 
gives I s =44GW/cm 2 , which is exactly the numerically 
found value (see Fig. [l]). Our numerical simulations con- 
firm this accuracy for all values of the domain length 
satisfying Ld ~ L c -C L and also for negative Ak. 

A first requirement for the QPM-induced separatrix to 
be relevant for switching purposes is that the intensity I s 



is low, which requires the ratio /3/7i to be small, i.e. the 
effective mismatch must be small (but nonzero) and the 
induced cubic nonlincarity strong. We will not discuss 
more appropriate materials than LiNb.Q 3 , but we note 
that the average model is very generaloEl, allowing 7! to 
be significantly increased by, e.g., modulating Ihe QPM 
gratingp or by a strong dc-value of the gratingp. 
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Nonlinear phase-shift of the FW versus inten- 
sity I for unseeded SHG in bulk LiNb0 3 . The QPM domain 
length is Ld=8.8fim. The crystal length is given in the figure. 
The vertical dotted line marks the separatrix 7 s =44GW/cm 2 . 
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tion (|7|) versus intensity. The minimum period between 
1=0 and the separatrix I=I S provides a good measure of 
the optimal crystal length observed in Fig. |2|, for which 
the 7r-shift due to the separatrix is the clearest and the 
flat phase versus intensity plateaus on either side are the 
broadest. From Fig. || we see how this optimal length in- 
creases as the domain length approaches the coherence 
length 8.9/zm, i.e. as exact effective phase-matching is 
approached. At the same time the separatrix intensity 
decreases. 

Thus it is desirable to work close to, but not exactly 
at, exact effective phase-matching in order to reduce the 
holding intensity in the switching process and keep a 
reasonable crystal length. Clearly the resolution in the 
photolithographic process could be an issue if the sepa- 
ratrix is to be measured and used. 

The research is supported by the Danish Technical Re- 
search Council through Talent Grant No. 5600-00-0355. 



So far the induced separatrix has not been observed in 
experiments. The reason could be: (i) The intensity was 
too low; (ii) the induced and inherent self-phase modu- 
lation (SPM) terms eliminate each other - conventional 
materials are self-focusing and have normal disper- 
sion and thus the inherent SPM coefficient is positive, 
while the induced SPM coefficient (—72) is negative; (hi) 
the crystal length was inappropriate. In Fig. |2| we illus- 
trate the effect of the crystal length. Clearly, the sepa- 
ratrix becomes more and more "hidden" in the overall 
variation due to the quadratic nonlinearity, and the flat 
plateaus on each side become narrower, as L is increased 
from the optimal value around 351/^m. When L is too 
short the separatrix is lost. 
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Fig. 3. Left: Period versus intensity I for unseeded SHG 
in bulk QPM LiNb0 3 with the QPM domain length given 
in microns at the curves. Right: Optimal crystal length and 
separatrix intensity 7 S versus domain length. 

In Fig. H we show the period 2K(m)/r of the solu- 
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